A Novel Time Integration Method for Solving A Large System of Non-Linear Algebraic Equations
نویسندگان
چکیده
Iterative algorithms for solving a nonlinear system of algebraic equations of the type: Fi(x j) = 0, i, j = 1, . . . ,n date back to the seminal work of Issac Newton. Nowadays a Newton-like algorithm is still the most popular one due to its easy numerical implementation. However, this type of algorithm is sensitive to the initial guess of the solution and is expensive in the computations of the Jacobian matrix ∂Fi/∂x j and its inverse at each iterative step. In a time-integration of a system of nonlinear Ordinary Differential Equations (ODEs) of the type Bi jẋ j + Fi = 0 where Bi j are nonlinear functions of x j , the methods which involve an inverse of the Jacobain matrix Bi j = ∂Fi/∂x j are called “Implicit”, while those that do not involve an inverse of ∂Fi/∂x j are called “Explicit”. In this paper a natural system of explicit ODEs is derived from the given system of nonlinear algebraic equations (NAEs), by introducing a fictitious time, such that it is a mathematically equivalent system in the n + 1-dimensional space as the original algebraic equations system is in the n-dimensional space. The iterative equations are obtained by applying numerical integrations on the resultant ODEs, which do not need the information of ∂Fi/∂x j and its inverse. The computational cost is thus greatly reduced. Numerical examples given confirm that this fictitious time integration method (FTIM) is highly efficient to find the true solutions with residual errors being much smaller. Also, the FTIM is used to study the attracting sets of fixed points, when multiple roots exist. 1 Department of Mechanical and Mechatronic Engineering, Taiwan Ocean University, Keelung, Taiwan. E-mail: [email protected] 2 Center for Aerospace Research & Education, University of California, Irvine Keyword: Nonlinear algebraic equations, Iterative method, Ordinary differential equations, Fictitious time integration method (FTIM)
منابع مشابه
Preconditioned Generalized Minimal Residual Method for Solving Fractional Advection-Diffusion Equation
Introduction Fractional differential equations (FDEs) have attracted much attention and have been widely used in the fields of finance, physics, image processing, and biology, etc. It is not always possible to find an analytical solution for such equations. The approximate solution or numerical scheme may be a good approach, particularly, the schemes in numerical linear algebra for solving ...
متن کاملUsing shifted Legendre scaling functions for solving fractional biochemical reaction problem
In this paper, biochemical reaction problem is given in the form of a system of non-linear differential equations involving Caputo fractional derivative. The aim is to suggest an instrumental scheme to approximate the solution of this problem. To achieve this goal, the fractional derivation terms are expanded as the elements of shifted Legendre scaling functions. Then, applying operational matr...
متن کاملThe Sine-Cosine Wavelet and Its Application in the Optimal Control of Nonlinear Systems with Constraint
In this paper, an optimal control of quadratic performance index with nonlinear constrained is presented. The sine-cosine wavelet operational matrix of integration and product matrix are introduced and applied to reduce nonlinear differential equations to the nonlinear algebraic equations. Then, the Newton-Raphson method is used for solving these sets of algebraic equations. To present ability ...
متن کاملAn Efficient Numerical Algorithm For Solving Linear Differential Equations of Arbitrary Order And Coefficients
Referring to one of the recent works of the authors, presented in~cite{differentialbpf}, for numerical solution of linear differential equations, an alternative scheme is proposed in this article to considerably improve the accuracy and efficiency. For this purpose, triangular functions as a set of orthogonal functions are used. By using a special representation of the vector forms of triangula...
متن کاملJacobi Operational Matrix Approach for Solving Systems of Linear and Nonlinear Integro-Differential Equations
This paper aims to construct a general formulation for the shifted Jacobi operational matrices of integration and product. The main aim is to generalize the Jacobi integral and product operational matrices to the solving system of Fredholm and Volterra integro--differential equations which appear in various fields of science such as physics and engineering. The Operational matr...
متن کاملHaar Matrix Equations for Solving Time-Variant Linear-Quadratic Optimal Control Problems
In this paper, Haar wavelets are performed for solving continuous time-variant linear-quadratic optimal control problems. Firstly, using necessary conditions for optimality, the problem is changed into a two-boundary value problem (TBVP). Next, Haar wavelets are applied for converting the TBVP, as a system of differential equations, in to a system of matrix algebraic equations...
متن کامل